Field theory analysis of S = 1 antiferromagnetic bond-alternating chains 

in the dimer phase 
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Dynamics of S — 1 antiferromagnetic bond-alternating chains in the dimer phase, in the vicinity 
of the critical point with the Haldane phase, is studied by a field theoretical method. This model 
is considered to represent the compound NXCgH^N^NC^ClCU (abbreviated as NTENP). We 
derive a sine-Gordon (SG) field theory as a low-energy effective model of this system, starting from 
a Tomonaga-Luttinger liquid at the critical point. Using the exact solution of the SG theory, we 
give a field theoretical picture of the low-energy excitation spectrum of NTENP. Results derived 
from our picture are in a good agreement with results of inelastic neutron scattering experiments on 
NTENP and numerical calculation of the dynamical structure factor. Furthermore, on the basis of 
the obtained theoretical picture, we predict that the sharp peaks correspond to a single elementary 
excitation are absent in the Raman scattering spectrum of NTENP in contrast to the inelastic 
neutron scattering spectrum. 



PACS numbers: 75.10.Pq, 03.70.+k, 75.40.Gb, 75.10.Jm 



I. INTRODUCTION 

One dimensional antiferromagnetic Heisenberg chains 
show various physical properties under strong quantum 
fluctuations. They have offered many fascinating phe- 
nomena as a stage of both theoretical and experimen- 
tal studies for many years. While an uniform antifer- 
romagnetic Heisenberg chain is gapless for half-integer 
spins, it has an exotic ground state with a finite gap 
for integer spins X This gapped phase is called Haldane 
phase. The antiferromagnetic Heisenberg chain is not 
exactly solvable except for 5 — 1/2. However, the Hal- 
dane phase could be understood in terms of the solvable 
Affleck-Kennedy-Lieb-Tasaki model^. Its exact ground- 
state can be constructed in terms of valence bonds (sin- 
glet pairs of constituent 5 = 1/2's). 

Bond-alternating antiferromagnetic Heisenberg chain 
defined by 



n = jJ2( S 2j-fS: 
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>2j • 02j+l 



), 



(1) 



where a represents the bond alternation ratio, is an inter- 
esting generalization of the Haldane gap problem. In this 
model, for spin quantum number S, the system shows 25 
successive quantum phase transitions in < a < oo^. 
This may be understood, in the AKLT picture, as fol- 
lows. In the limit of a = 0, all the valence bonds are 
placed on the "odd" links between sites 2j—l, 2j. In the 
opposite limit a — oo, all the valence bonds are placed 
on the "even" links between sites 2j, 2j + 1. They corre- 
spond to completely dimerized states. By changing the 
number of valence bonds on odd and even links, we can 
construct 2S+ 1 AKLT-type states with varying degree of 
dimerization. Each of these AKLT-type states represents 



the 25+1 gapped phases, separated by the 2S quantum 
phase transitions. Each of the transitions may be re- 
garded as a rearrangement of valence bonds, in which one 
valence bond is transferred between neighboring links. 

In the case of an integer spin 5, 25+1 gapped phases 
can be classified into two categories: one with an even 
number of valence bonds on every link, and the other 
with an odd number of valence bonds. The latter is a 
(symmetry-protected) topological phase; in the presence 
of global Z2 x Z2 symmetry of 7r-rotation about x, y and 
z axes, it is characterized by spontaneous breaking of a 
hidden Z 2 x Z 2 symmetry^, which can be detected by 
the string order parameter^ More generally, the topo- 
logical phase is also protected by either the time-reversal 
symmetry or the lattice inversion symmetry about a link, 
and is characterized by exact two-fold degeneracy of en- 
tire entanglement spectrumi&iS 

In this paper, we focus on the simple case of 5 = 1. 
Here, the 25 + 1 = 3 gapped phases consist of two dimer 
phases (at large and small values of a) and the 5 = 1 Hal- 
dane phase in a range of a including a = 1 . As mentioned 
above, being an exotic "topological phase" , the Haldane 
state has been studied vigorously. In contrast, the dimer 
phase may be regarded as a trivial phase as far as the 
ground-state properties are concerned, and has been rela- 
tively less studied. In fact, the system is exactly solvable 
in the dimer limit a = 0, where the system is reduced 
to a two-body problem. However, dynamical properties, 
which are related to excited states, are not necessarily 
trivial even in the dimer phase. While quantum phases 
are usually classified with respect to order parameters in 
the ground state, they are not sufficient for elucidation 
of the excited states. Thus the dimer phase, although 
without any nontrivial order, has a possibility of exhibit- 
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ing a rich structure in the excitations. In addition, it 
would be also interesting to compare dynamical prop- 
erties between the dimer phase and the Haldane phase. 
Therefore, in this paper, we attempt to study the dy- 
namical properties in the S = 1 dimer phase, especially 
in the neighborhood of the quantum critical point. 

Theoretical study of Haldane gap has been, without 
doubt, stimulated by discovery of several model ma- 
terials, such as Ni(C 2 H 8 N 8 ) 2 (N02)C104 (abbreviated 
to NENP) and Ni(C 5 H 14 N2)2N3(PF 6 ) (abbreviated to 
NDMAP). Fortunately, materials representing the 5=1 
bond-alternating antiferromagnetic chain ((T|), such as 
Ni(C 9 H 2 4N 4 )(N02)C104 (abbreviated to NTENP)Ii and 
[Ni(333-tet)(^-N 3 )„](C104)„ were also discovered, and 
the bond-alternation ratio a was identified in these mate- 
rials. In the latter system, quite remarkably, a turns out 
to be exactly (within experimental accuracy) at the crit- 
ical point separating the Haldane and dimer phases. 12 In 
NTENP, on the other hand, a is slightly smaller than the 
critical valued That is, the system belongs to the dimer 
phase, but in vicinity of the quantum critical point. Thus 
NTENP offers an ideal playground to study nontrivial 
dynamics in the dimer phase. 

Dynamical properties of NTENP are studied experi- 
mentally by inelastic neutron scattering (INS)$^~— and 
arc numerically studied by using a continued-fraction 
method based on the Lanczos algorithm 17 . However, in 
order to better understand the dynamics in NTENP and 
also in related materials in a unified way, it would be 
desirable to develop a coherent theoretical picture. In 
this paper, we present a field theoretical approach to the 
problem, which we hope would serve the purpose. 

At a quantum critical point, physical quantities are 
often constrained by the conformal invariance, and the 
system is described by a conformal field theory (CFT). 
Near the critical point, the system could be understood 
by a CFT with a small (but relevant) perturbation. In 
this paper, we follow this strategy, making use of the fea- 
ture of NTENP that the bond-alternating ratio is in the 
vicinity of the critical point. Based on the obtained ef- 
fective theory, we elucidate the dynamical properties of 
NTENP and discuss preceding experimental and numeri- 
cal results. Furthermore, we make predictions on Raman 
Scattering experiment which has not yet been carried 
out. 



II. MODEL AND METHOD 

In fact, NTENP has a sizable spin anisotropy; it is bet- 
ter described, instead of eq. ([1} , by the following Hamil- 
toniani^ 

H = J£(S 2 ;-i ■ S 2j +aS 2j ■ S 2J+1 ) + DJ2(.S*) 2 , (2) 

3 3 

where J > 0, a and D represent bond alternation and 
uniaxial anisotropy, respectively. For NTENP, these pa- 
rameters are D ~ 0.25J and a ~ 0.45. The spin 



anisotropy axis (which is determined by the chain axis) 
is taken to be z-direction. In this paper, we ignore the 
in-plane anisotropy EJ2j(Sj) 2 — (S^) 2 and the exchange 
anisotropy, which are expected to be smaller than D in 
NTENP. The lattice constant between neighboring two 
spins is set to unity. 

The critical point a c , where energy gap is closed, is 
about 0.6 at D = 0. It is a function of D which is numer- 
ically evaluated in Ref. [H. According to them, a c is al- 
most constant for < D < 0.3 J. In NTENP, D ~ 0.25J 
and thus we can assume that a c ~ 0.6. The actual bond- 
alternating ratio of NTENP is a ~ 0.45 < a c , and the 
system is certainly in the dimer phase. However, it can be 
still regarded as a vicinity of the quantum critical point. 

At the critical point a = a c , the system is mapped 
to a TL liquid, 19 in the low-energy limit. 13 It is nothing 
but the theory of free boson field </> with the Lagrangian 
density, 

£=\{d»cp) 2 - (3) 

Let R be the compactification radius of boson field, 

<f>(x,t) + 2irR~ (f>(x,t). (4) 

The compactification radius determines critical expo- 
nents of the TL liquid. This may be understood as a con- 
sequence of "quantization" of vertex operators required 
by the compactification. For example, the vertex oper- 
ators cos ('yip) must have the coefficient j = n/R with 
an integer n, in order to be single- valued under the com- 
pactification (|4]). The scaling dimension of the vertex 
operator cos (n<f>/ R) is known^ as 



The operator is relevant (irrelevant) in the 
Renormalization-Group (RG) sense, if the scaling 
dimension is smaller (greater) than 2. Thus, among 
the vertex operators of the above form, n — 1 is the 
most relevant. As we will discuss later, in the present 
application R is somewhat smaller than the SU(2) 
symmetric value l/y2w, implying n = 1 is the only 
relevant operator in the family cos (n<j>/R). 

When a y^z a c , the system acquires a gap and it is no 
longer described by the TL liquid (|3|). However, if \a— a c \ 
is small, as it is the case in NTENP, the low-energy ef- 
fective theory would be given by the TL liquid with a 
perturbation. We postulate that the leading perturba- 
tion is the most relevant operator permitted under the 
symmetry of the Hamiltonian ([2]) : cos (<p/R). Hence the 
low-energy effective theory for small \a — a c \ should be 
following sine-Gordon (SG) field theory, 

£=^) 2 + Ccos(|), (6) 

where C is the constant that is proportional to a — a c 
and R is the compactification radius of boson field at the 
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critical point. When a = a c , C equals zero and Eq. (J6j) is 
reduced to Eq. ([3]). It should be noted that we use R at 
the critical point a c in this construction (j6j for a 7^ a c . 
The compactification radius R is well defined only at the 
critical point and the precise determination of R is subtle 
for a 7^ a c . Nevertheless, R obtained at a = a c should 
be a reasonable approximation as long as \a— a c \ is small. 
We also note that R at the critical point a c is a function 
of D. 

The SG field theory is integrable and its exact solu- 
tion is known j 20 i 21 In this study, we analyze the low- 
energy spectrum of NTENP using the exact solution of 
the SG field theory, where the original spin Hamiltonian 
is mapped. We thus refer to this exact solution in the 
next section. 



III. EXACT SOLUTION OF THE 
SINE-GORDON FIELD THEORY 

The SG field theory has the Lorentz invariance, hence 
the energy of the elementary excitations obey the rela- 
tivistic dispersion relation 



\J k 2 v 2 



(7) 



Here, v s is the spin-wave velocity which plays the role of 
the speed of light, k is the momentum, and m is the mass 
of the excitation (in the unit of energy) , which represents 
the lowest creation energy of each elementary excitations. 

The elementary excitations of the SG field theory © 
consist of a soliton, an antisoliton and several breathers. 
The soliton and the antisoliton have the same mass Ms- 
Breathers are bound states of a soliton and an antisoliton. 
The exact solutio n 20 i 21 of the SG field theory implies the 
number of different kinds of breathers as [l/£], where 



- = 8irR 2 - 1, 

and [x] denotes the integer part of x. 
exact mass ratio was also derived as 



(8) 



Remarkably, the 



= 2 sin I — 1 

M s y 2 ' 



(n < 1/0, 



(9) 



where M n and Ms are respectively the n-th breather 
mass and the soliton mass. Here we note that the 
breather mass M n does not exceed twice the soliton mass 
2M5. Physically, this can be interpreted as follows. If the 
mass of a breather were larger than 2Ms, the breather 
would decay into free soliton and antisoliton, and would 
not constitute a stable elementary excitation. It should 
be emphasized that the mass ratio (j9|) and the number of 
different breathers depend only on the boson field com- 
pactification radius R, which is a function of D. 

The above picture of elementary excitations implies the 
spectrum of the excited states schematically described in 
Fig. [1] Each elementary excitation obeys the relativistic 
dispersion ([7]) with a species-dependent mass m and the 
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FIG. 1. Schematic view of the excitation spectrum of the 
SG field theory. |GS), \S), \S), \Bi) and \B 2 ) represent the 
ground state, the soliton state, the antisoliton state, the first 
breather state and the second breather state, respectively. 
The soliton level and the antisoliton level is degenerated. The 
number of different breathers is [l/£], where [x] denotes the 
integer part of x and l/£ = — 1. In this schematic view, 
the case that the theory has two breathers is represented. 
Owing to the Lorentz invariance of the theory, the excitation 
energies obey the relativistic dispersion relation (FT}. 



common spin- wave velocity v s . Above the energy 2Ms, 
there is a continuum of excited states, corresponding to 
scattering states of two or more elementary excitations. 

The elementary excitations of NTENP is understood 
as the corresponding elementary excitations of the SG 
field theory. The boson field compactification radius R 
of NTENP determines the mass ratio between a soliton 
and breathers and the number of different breathers cor- 
respond to NTENP. 

In passing, we note that, for an isotropic system, 
the SU(2) symmetry requires R = 1/v2tt- The exact 
solution of the SG theory implies that there arc two 
breathers, and the mass of the 1st breather is degenerate 
with soliton/antisoliton, forming a SU(2) triplet. Fur- 
thermore, the mass of the second breather is y/3 times 
that of the triplet. However, this does not quite agree 
with the numerical results. This is due to the marginal 
perturbation to the SG theory, which generally exists 
in the SU(2) symmetric system. Only with an appro- 
priate strength of next-nearest neighbor interaction, the 
marginal operator can be eliminated and the above pre- 
diction of the SG theory is realize d 22 ' 23 . In this pa- 
per, however, we consider a system with a substantial 
anisotropy, where the marginal operator is absent. Thus 
the simple SG theory ^ is expected to be a good low- 
energy effective theory of the system. 
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IV. STRUCTURE OF ELEMENTARY 
EXCITATIONS 

A. Mass ratio among the elementary excitations 

We are now in a position to apply the relation © to the 
system of NTENP. Here we note again that the critical 
point a c and the compactification radius R at this point 
are a function of the uniaxial anisotropy D. At D = 0, 
R is determined as 2ttR 2 = 1 by the SU(2) symmetry of 
the system. The variation of R at several values of D is 
already studied numerically by using CFT and the level 
spectroscopy method, by Chen et al*i£ The parameter K, 
5 and D in Ref. HI correspond to l/{2nR 2 ), (l-a)/(l + 
a) and 2D/(1 + a) respectively, in the notation of the 
present paper. Thus, for NTENP with D ~ 0.25J, we 
find l/(2nR 2 ) ~ 1.7 that implies l/£ - 1.2. It follows 
that this effective SG field theory for NTENP has the 
first breather only. The breather/soliton mass ratio is 



Mi 

M s 



(10) 



Accordingly, the energy spectrum of NTENP based on 
the SG field theory picture appears to be the following: 
the lowest excited state is a degenerate doublet corre- 
sponds to the soliton and the antisoliton, the next lowest 
excited state is a single mode corresponds to the breather, 
and an excitation continuum starts from the energy twice 
the soliton mass. This spectrum is schematically de- 
scribed in Fig. 01 



B. Dynamical structure factor 

The dynamical structure factor (DSF) is defined as 
the Fourier transform of the correlation function. It is 
an important quantity, as the differential cross section 
in an INS is proportional to the DSF. The fi (= x, y, z) 
component of the DSF at zero temperature is 



S»»(q,u)= I dtJ2{i>o\S?(t)S%(0)\i> )e 



iu)t—iq(j — k) 



3G 



= J2 S ^- ( E » - E °)) KWlVo}| 2 , (ii) 

n 

where |V>o) and \i/) n ) are the ground state and the excited 
state of Hamiltonian ([2} with the eigenvalue Eq and E n , 
respectively, and Sjf = -j== Yfj=i S^e~ tqj , where 2N is 
the total number of spins. Here it should be noted that, 
in the systems with a bond alternation such as NTENP, 
the width of the Brillouin zone is reduced to half of that 
of the uniform chain. Namely, the conserved crystal mo- 
mentum can be defined as 

q = q mod 7T, (12) 

which is defined on the reduced Brillouin zone 

-7r/2<q<7r/2 (13) 



In particular, q = and q — tt are both identified with 
the center of the reduced Brillouin zone q = 0. However, 
it should be noted that, the DSF is not equivalent for q 
and q + tt. This could be easily understood by consid- 
ering the limit of zero dimerization, where q and q + tt 
are certainly distinguishable. The operator and S^ +7T 
are different operators and thus give different DSFs. In 
the presence of dimerization, they are governed by the 
same selection rule. The matrix elements, however, are 
generally different. 

The summation in Eq. (1111) is taken over all inter- 
mediate states, but many of the matrix elements vanish 
due to selection rules. While excited states \ip n ) contain 
one or more elementary excitations, only the states which 
have non- vanishing overlap with S£\?po) contribute to the 
DSF. Thus, in order to evaluate each component of the 
DSF, it is important to know what kind of elementary 
excitations are created by applying on the ground 
state \ipo)- F° r example, if S^\tpo) contains a single soli- 
ton state, there must be a contribution proportional to 
S(u> ~ ^q 2 + Mf), to S^{q,uj). 



C. Dimer limit 

In order to clarify which elementary excitation is cre- 
ated by application of each spin operator, we need to 
know the correspondence between the operators in the 
SG field theory and those in the original spin system. In 
principle, this would follow from a microscopic derivation 
of the effective SG field theory from the spin model ([2]). 
However, the microscopic derivation is rather compli- 
cated and the correspondence cannot be easily estab- 
lished. A possible derivation starts from the bosonization 
of two S = 1/2 chains and then introduces a strong ferro- 
magnetic coupling between them, to form aS = l chain 
effectively^. Instead of pursuing this direction, we shall 
develop a simple ansatz for the correspondence, based on 
the dimer limit a = 0. 

In the dimer limit, the problem is reduced to a two-spin 
problem, and thus can be exactly solved. Two S = l's 
coupled with the Heisenberg antiferromagnetic exchange, 
H = JSi ■ S2 has the singlet groundstate with energy 
—2 J, triplet excited states with energy — J, quintet ex- 
cited states with energy +J. Now let us introduce the 
uniaxial anisotropy D to the dimer, and consider the 
Hamiltonian 

7i = JS 1 -S 2 +D({S Z 1 ) 2 + (S Z 2 ) 2 ). (14) 

The singlet groundstate is perturbed by the anisotropy, 
and given by 



\s) =CJ\ U) + I It) + 



A - VA 2 + 8 



|0G) , (15) 



with the groundstate energy 



A-VA^T. 



■J, 



(16) 
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where A = 2D — 1 and C s is the normalization constant. 

The uniaxial anisotropy splits the triplet and quintet 
excited states. Each state can be labelled by the to- 
tal magnetization S z — S* + 3%, which is still a good 
quantum number. Since the uniaxial anisotropy does not 
break the 7r-rotation symmetry under rotation about x- 
axis. Thus, the triplet states are split into a singlet with 
S z = 0, which is denoted by |t°), and a doublet 1^) with 
S z = ±1. The doublet states are given by 

l* + ) = ^(lt0)-|0t)), (17) 

\t-) = ±(\i0)-\0i)), (18) 

with the energy D — J, while the singlet is given by 

|f°> = ^|(in>-ut», as) 

with the energy 2D — J. For D > 0, which is the case 
for NTENP, the doublet states j^) are lower in energy 
than the singlet \t°). 

In a similar way, the quintet states are split into a sin- 
glet with S z = and two doublets with S z = ±1 and 
S z = ±2. However, as we will argue below, once we 
move away from the dimer limit, they decay into excita- 
tions corresponding to |t 0,± ). Thus, for the purpose of 
classification of elementary excitations, it is sufficient to 
consider the (split) triplet |t 0,± ). 

Given these states in single dimers, the states in the 
entire system is given as follows: The ground state is 
the state that all dimers form \s). The lowest excited 
state is the state that \t + ) or \t~) is excited in one of the 
dimers, while the other dimers remain in \s). Similarly 
in the next lowest excited state, \t°) is excited in one 
of the dimers and the other dimers form \s). Here, we 
denote these states |0), |one t + ), |one t~), and |one t°), 
respectively. In fact, in a system which consists of N 
dimers (27V sites), there are N degenerate states for each 
class of |one i^), where £ = 0, ±, corresponding to the 
location of the excited dimer. Let us define lone 
where £ = 0, ±, as the state with the j-th dimer (of sites 
2j — 1 and 2j) as the only excited dimer. We can also 
define their Fourier transform 

|one ti, q) = -±=J2 ^{one t<,j), (20) 
j 

where C — 0, ±, labeled by the crystal momentum q de- 
fined in the reduced Brillouin zone (|13l) . 

In this paper, we are primarily interested in the DSF 
around the antiferromagnetic wavevector, q ~ ir. Here, 
the matrix elements of spin operators among these states 
are found to be: 

(onet+,q\S+\0) ^0, 
(onet-,q\S-\0)^0, (21) 
(one t°,q\S z \0) ^0, 



generically if q = q mod 7r, and 

the other matrix elements of S£ 
between |0) and |one t^), or |one t°) = 0. 

Now let us consider increasing a from zero. With a 
nonzero a, an excited state in a single dimer can "hop" 
to the neighboring sites, giving rise to a dispersion. It 
should be noted that, in the dimerized limit, the momen- 
tum basis (f20| as well as the localized basis |one 
represents N degenerate excited states. However, when 
the interdimer interaction is introduced with a > 0, mo- 
mentum basis ([20| , but not the localized basis gives a set 
of approximate eigenstates if a is sufficiently small. Now 
the energy of the excited state |onei'», q) depends on the 
momentum q; the dependence is nothing but the disper- 
sion relation. Since the structure of the triplet split into 
the doublet and the singlet by the anisotropy is very sim- 
ilar between the single dimer and the SG field theory, it 
would be natural to assume that the structure of the ex- 
citation does not qualitatively change in < a < a c . Un- 
der this assumption, the states in the dimer limit a = 0, 
|0), |one t + ), |one t~) and |one t°) are considered to be 
smoothly connected to, at a ~ a c , the ground state |GS), 
and the single soliton state \S, q), the single antisoliton 
state \S, q), the single breather state \B\,q), each with 
the momentum q. If this is the case, the following selec- 
tion rules should hold: 

(S,q\S+\GS)^0, 

(5,«|5-|GS)^0, (22) 
(B 1 ,q\S z q \GS)^0, 

generically if q = q mod 7r, and 

the other matrix elements of S£ 

between |GS) and \S,q), \S,q), or \B u q) = 

Assuming this, S^, 5~ and S z respectively act as the 
creation operators of the soliton, the antisoliton and the 
breather. Thus the soliton and the antisoliton contribute 
to the x, y components of the DSF and the breather does 
to the z component. 

The assumption above is natural because we are in- 
terested in a system in the dimer phase, and there is no 
phase transition between the dimer limit a = 0. How- 
ever, properties of excited states are not necessarily the 
same, even if they are in the same phase and the ground- 
states are adiabatically connected. We confirmed the va- 
lidity of our assumption and that the matrix elements in 
Eq. (f2"2"]l does not vanish, by the numerical exact diag- 
onalization calculation using the Lanczos method. We 
show our numerical results in Fig. [2] and Fig. [3l In what 
follows, we direct our attention to the states with q = 
which correspond to the bottom of the each branch in 
the energy spectrum, in order to compare with the pre- 
ceding results i 14 i 15 i 17 As we will clarify later, it suffices 
to obtain the matrix element of the local spin operator 
between the groundstate and the excited states. 
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FIG. 2. (a) The a-dependence of the energy of the ground 
state and the lowest excited state in subspaces that the total 
magnetization m — . Sj is ±1, and the momentum q = 0. 
(N = 16) The case that a = 0.45 corresponds to NTENP. 
(b) The excitation energies extrapolated to N — >• do. Inset: 
Extrapolation of the excitation energies at a = 0.45 to TV — > 



Fig. [U(a) shows the energy of the ground state and the 
lowest excited state in subspaces that the total magne- 
tization m = Y]j Sj is ±1,0 and the momentum q = 0. 
This figure is obtained by the exact diagonalization of 
the 16 site spin system. Fig. [2] (b) shows the exci- 
tation energies extrapolated from the calculated values 
in the 12, 14 and 16 site spin systems to the infinite 
size. These results show that the energy of the elemen- 
tary excitations changes adiabatically under the varia- 
tion of a from to 0.45 (~ a c ), which is the value of 
NTENP. Thus our assumption that the excitation struc- 
ture evolves smoothly under the variation of a, is con- 
firmed. In what follows, we represent the states at a that 
are adiabatically connected from \<p) at a = 0, as \4>)a- 
The states |one t + ) a ~a C ) |one t~) ar ^ ac and |one t z ) a ^ ac 
correspond to the soliton | S) , the antisoliton | S) and the 
breather \B\), respectively. When a is finite, the elemen- 
tary excitations have a dispersion. Therefore we label 
them also with the momentum q. 

Fig. [3] shows the absolute value of the matrix ele- 



,(one t+\Sf\o) a , „<one t~\SJ\0) a 
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FIG. 3. The a-dependence of the absolute value of the ma- 
trix elements a (one t + ,q = 0\Sf \0} a , a (one t~,q = 0|S~|0) C 
+o - 



and a (one t",q = 0\S-\0) a . 
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16) They correspond to 



(S|Sr|GS), (S|S.-|GS) and (B^S^GS) in the vicinity of the 



critical point a — a c , respectively. Inset: 
the matrix elements at a = 0.45 to N — > cx 



Extrapolation of 



,q = 0\ST\0) a 

one t°,a — 0\Sj\0) a in the variation of a. This 



ments Q (one t + ,q = 0|Sy~|0) a , Q (one t 
and Q (one t°,q - , , 

figure is also obtained by the 16 site spin system ex- 
act diagonalization. From these calculated values, the 
matrix elements of S£ corresponding to them are eval- 
uated as follows: As we will discuss later in Sec. I VII 
the ground state |0) is link-parity even, and the states 
|one t + ,q = 0) Q , |one t~,q = 0) a and |one t z ,q — 0) a are 
link-parity odd: 



P/|onc t c ,q = 



Pi\0) a = |0) a , 

Q) a = -|one t c ,q = 0) Q , 



where C = 0, ±, and the link parity operator Pi represents 
a reflection with respect to a link. It satisfies Pf = I, 
where / is the identity operator. From these properties, 
it follows that 

Q (one t c ,q = 0\Sf +1 \0) a = a (one ,q = 0\P t S+ Pi\0) a 

= - a (one t<,q = 0\S+\0) a , 



where ^ = 0,±. Accordingly, , 
does not vanish, while „ (one , 



(23) 

(one *C,g = 0|S+ = J0) Q 
,g = 0|S'+ = o|0) Q =0. The 
adiabatic continuity of the excited states then implies the 
same selection rule for the single soliton state \S, q = 0), 
single antisoliton state \S,q — 0), and single breather 
state \Bi, q = 0). Thus we have confirmed the validity of 
our ansatz. 



D. 



Excitation structure of NTENP based on the 
SG field theory 



Combining the results in Sees. IIV Al and [TV Bl we ob- 
tain the following picture: NTENP has three elementary 
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E" 



2M S 
1.93M S 



M S 



continuum 

\Bi) 

\s)AS) 

|GS> 



citations, then only those of Sjj =Tr contribute. For exam- 
ple, the contribution of the soliton with q — occurs in 
S xx (tt, u)s)- Taking these into consideration, the inten- 
sity ratio are evaluated by 







FIG. 4. Schematic view of the excitation spectrum of NTENP 
based on the SG field theory. The doubly degenerated isolated 
mode correspond to the soliton and the antisoliton is present 
and its lowest creation energy is Ms- The system has only the 
first breather and its lowest creation energy is 1.93Afs. The 
excitation continuum starts from the 2Ms at the momentum 
g = 0. 



excitations that correspond to the soliton, the antisoli- 
ton and the first breather in the SG theory. For these 
elementary excitations, the soliton mass and the antisoli- 
ton mass are degenerate, and the first breather mass is 
1.93 times the soliton mass. (The mass degeneracy of 
soliton/antisoliton is a common nature of the SG field 
theories, that is, irrespective of the value of radius R. 
On the other hand, the breather mass ratio depends on 
the value of the radius R. For NTENP, 1/(2ttR 2 ) ~ 1.7.) 
Thus, it is concluded that the doubly degenerated iso- 
lated mode is at Ms and the single isolated mode is at 
1.93Afs in the energy spectrum of NTENP. The excita- 
tion continuum starts from 2M$- This is schematically 
described in Fig. |4j 

The soliton, the antisoliton and the breather are re- 
spectively created by operating S+, S~ and £* on the 
ground state. Conversely, when each spin operator cre- 
ates a single elementary excitation, its species is uniquely 
determined as soliton, antisoliton, and breather. Hence 
the single-particle peak due to the soliton/antisoliton 
contribute only to S xx and S yv , and that due to the 
breather does to only S zz . The intensity ratio between 
the lowest excited states in S xx and S zz are evaluated 
from the result in Fig. [3] Since the bottom of each ele- 
mentary excitation branch in the energy spectrum q = 
corresponds to q = and q = ir, the intensity from the 
lowest excited state in <5> pp is the sum of S^fOjO;) and 
l S At/i (7r, lo) (where oj is the energy transfer of the lowest 
excited state). However, the matrix elements of S^ =0 
vanishes due to the relation (|23p for each elementary ex- 



KBr^Ol^JGS)! 5 



r|(S,<? = 0|S+=.|GS)| 2 + \\(S,q = 0|S^|GS)| 2 



0.22. 

(24) 



Here, we use the calculated values .45( one t + ,q = 
0|5'+ =7r |0)o.45, o.45 (one t~,q = 0|S , g=7r |0) . 4 5 and 
o.45(one t z ,q — 0|S'| =w |0)o.45 in Fig. [3] for the soliton, 
the antisoliton and the breather, respectively. 



V. COMPARISON WITH PRECEDING 
RESULTS 

According to Refs. [HI and [TU, two main peaks are 
present in the INS spectrum of NTENP. The first peak 
occurs from the fluctuation of spin x, y components and 
the second from that of z component. In an inelastic 
scan at the center of the Brillouin zone, the energy ra- 
tio between the peak from spin x, y component and that 
from z component is about 1.91/1.07 ~ 1.79, and the 
intensity of the second peak is about 0.2 times that of 
the first peak. Furthermore, the DSF in the model ([2]) 
under a transverse magnetic field H is numerically stud- 
ied using the continued-fraction method^. Their result 
at H = can be compared with the present study. The 
ratio between the energies correspond to the first peak 
in S xx (ir, oj) and S zz (tt,uj) is ~ 1.91 and the intensity 
of the first peak of S zz (tt,lu) is about 0.22 times that 
of S xx (ir,u). We also note that, the ratio of the lowest 
gap to m = ±1 and m = excitations in our exact di- 
agonalization study, shown in Fig. (5J is consistent with 
the predicted mass ratio 1.93. This gives an additional 
support to our field-theory analysis. 

The lowest-energy peak in the INS spectrum in Refs. [13 
and [H, which corresponds^ to the lowest excited state in 
S xx (ir, uj), is the contribution of the soliton/antisoliton in 
our picture. Likewise, the second peak in the INS spec- 
trum, which corresponds to the lowest excited state in 
S zz (n,uj), is the contribution of the breather. Our re- 
sults that the mass ratio between soliton and breather 
Mb /Ms ~ 1.93 and that the creation operator of soli- 
ton, antisoliton and breather is respectively S+, S~ and 
S*, agree very well with these results. Our result (|2~4")l 
on the intensity ratio also shows a good agreement. Not 
only the peak positions but also the polarization of the 
DSF is consistent. Our field-theory approach thus pro- 
vides a coherent theoretical picture describing the pre- 
ceding experimental and numerical results, at least at 
zero magnetic field. A slight disagreement with experi- 
ments may be attributed to the following factors. First, 
we neglected the in-plane anisotropy E, which must be 
present in NTENP with a triclinic crystal structure. Fur- 
thermore, the low-energy asymptotic description based 
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on the SG field theory is not exact for NTENP due to its 
finite energy gap. 

More recently, ESR study of NTENP was reported 
in Ref. 25. We will make a brief comment on ESR in 
Sec. Em 



VI. RAMAN SCATTERING IN NTENP 

As a new application of the present approach, in this 
section, we discuss the Raman scattering (RS) spectrum 
of NTENP and predict its peak positions. Following our 
discussion in Sec. lIVBl we first identify the non- vanishing 
matrix elements in the dimer limit a = 0, and then ex- 
tend the results to a ^ using the continuity of the 
excitation structure from a = to a < a c . 

According to the theory of exchange-scattering in mag- 
netic compounds, the spectral function J(w) for the scat- 
tered light at zero temperature is given by — 

J(w) = 5^J(w - {E n - E ))\(i; n \HB.\ip )\ 2 , (25) 

n 

where \ijjo) and |^ n ) are the ground state and the excited 
state of the Hamiltonian ©, whose energy eigenvalues 
are Eq and E n , and T~Lr is the effective interaction for 
the exchange-scattering of light, 

Hr = y^X e i ' di : i+i){e s ■ d i: i +1 )(Si ■ S i+1 ), (26) 

i 

for the present one dimensional system @. Here, e,, 
e s are unit polarization vectors of incident and scattered 
light, and d i i+ i is a unit vector connecting the spin sites 
i and i+ 1. In Eq. (|26[) . di j+i is constant in the present 
system @; once the experimental apparatus is set up, 
and e s are fixed. (For a recent discussion of related prob- 
lems, see Ref. [27h Thus, we just need to know matrix 
elements of the operator . Si ■ Stj+i in order to calculate 
RS spectrum. This effective Raman Hamiltonian (|21)1) . 
taking the form of an inner products of spin operators, 
possesses several symmetries. Consequently, there are 
following conserved quantities: the total magnetization 
m = Y) j Sj, the momentum and the link parity of the 
states. Link parity transformation is a reflection with 
respect to a link. In the systems with bond alternation 
such as NTENP, the link parity is conserved but the site 
parity (reflection with respect to a site) is not conserved. 
The link parity operator Pi satisfies _P ; 2 = /, where / 
is the identity operator. As mentioned above, the se- 
lection rule for RS is quite different from that for INS 
represented in Eq. (fTTj). This difference leads to distinct 
peak positions between the INS and the RS spectrum of 
NTENP. 

In what follows, we discuss the overlaps between \ip n ) 
and |^o) to detect the excited states that contribute to 
the RS spectrum. For this purpose, here we collect the 
conserved quantities in the present system ^ again: the 



total magnetization m = ■ , the momentum q asso- 
ciated with the invariance under the translation by two 
sites, and the link parity. 

First we analyze the a = case. The ground state |0) is 
the eigenstate of m with the eigenvalue 0, while the lowest 
excited states |one t + ) and |one t~) are that with the 
eigenvalue 1 and —1, respectively. Since T-Lr conserves 
the total magnetization m, these excited states do not 
have overlap with T-Lr\Q), and do not contribute to the 
RS spectrum. Thus in what follows, we only discuss the 
subspace with m = in order to detect the states having 
non- vanishing overlap with Hr\0). The three low-energy 
excited states in this subspace, applied the same notation 
as the preceding one, are |one f°), |one t + ,one t~), and 
| two t°), from the lowest in energy, in the dimer limit 
a = 0. They correspond to the single breather state, the 
scattering state of a soliton and an antisoliton, and the 
scattering state of two breathers, respectively at a ~ a c . 

The link parity of the ground state |0) and these states 
are calculated as follows: For a given state \4>), we denote 
its mapping by the link parity as \(j>) = Pi\4>)- As a 
preparation, we first discuss the link parity of the states 
in a single dimer. The state \s) is link-parity even: 

\s) = p t \ s ) = c s Pt (| n> + 1 it) + - — irrfl 00 » 

A — V A 6 + 8 

= g.(Ut) + lW + A _^ a + 8 |oo)) 

= |*) s (27) 

where A = 2D — 1 and C s is the normalization constant. 
In a similar calculation, we find that the states \t+), 
and \t~) are link-parity odd: 

\&) = P l \t°) = -\t }, (28) 
\t+)=P l \t + ) = -\t + ), (29) 
|t->=fl|t-> = -|t->- (30) 

Based on these, we next discuss the states in the entire 
system. We note that Pi changes the position of the 
dimers in addition to the action on each dimer described 
above. The link parity of |0) is even: 

P l \0)=P l (\s) ••») 
= \s)---\s) 

= | S )...| S ) = |0). (31) 

Next we discuss the link parity of excited states in the 
subspace with m — 0. A state with a fixed position of 
the excited dimer is not an eigenstate of Pi nor of the 
momentum. For example, the state |one t°) with the 
first dimer excited is transformed under Pi as 

P l (\t°)\s)\s)---\s))=-\s)\s)---\s)\t°). 

On the other hand, momentum eigenstates can be con- 
structed as a linear combination of 

|t»| S )...| S ), | S )|t»...|s>> 
| S )| S )|i°)...| S ), |*>|*>».|*>l*°>- 
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Since the ground state |0) has zero momentum and %r 
is translationally invariant, only the zero-momentum ex- 
cited states contribute to the RS spectrum. We thus only 
have to examine the zero-momentum combination 

|onei°,g = 0) oc |t°)|s)|s) • • • \a) 

+ | a >|t»..» + ... 
+ \s)\ S )---\s)\t ). 



It turns out that this is link-parity odd: 



Pi |one t°,q 



0) 



|one t°, q ■■ 



0). 



(32) 



Since the ground state |0) is link-parity even, |one t°, q = 
0) appears not to have the overlap with Hr\0). Similarly, 
|one t + ,q = 0) and |one t~,q = 0) are link-parity odd, 
and do not have the overlap. 

In the above analysis, we considered the dimerized 
limit a — >■ 0. However, the Hamiltonian respects the 
link parity Pi for any a. Thus, the adiabatic continuity 
of the excited states, which was demonstrated numeri- 
cally in Sec. IIV CI implies that the corresponding states 
belong to the same eigenvalue of Pi for any a < a c . It 
follows that, none of the "single elementary excitation" 
states | one t**) a (£ = 0, ±) (corresponding to a single soli- 
ton, a single antisoliton and a single breather for a ~ a c ) 
contribute to the RS spectrum for a < a c . 

Thus only the states in the excitation continuum 
have the possibility of having non- vanishing overlap with 
Hr\0). In particular, the matrix element with the "two 
elementary excitations" state (one t + ,one i~|%i?,|0) does 
not vanish in the dimer limit. This transition is not for- 
bidden by the selection rule due to the link parity P;. 
Thus it is natural to expect that it remains non-zero 
for a < a c . Indeed, we examined the matrix element 
a (one i + ,one t~\HR\0) a numerically for a < a c . The 
result, obtained after extrapolation to infinite size from 
the systems with length 12, 14 and 16, is shown in Fig. [5] 
This confirms that the matrix element does not vanish 
within the dimer phase, as expected. The structure of 
excitations is continuous up to the critical point a = a c , 
also in this regard. 

Thus the "two elementary excitations" state 
| one t + ,one t~) a should contribute to the RS spectrum 
of NENP. In the vicinity of the critical point a ~ a c , 
this state may be regarded as the scattering state of a 
soliton and an antisoliton. 

From these results, we predict that the sharp peaks 
corresponding to creation of single elementary excitation 
are absent in the RS spectrum of NTENP. The RS spec- 
trum consists only of excitation continuum, and is quali- 
tatively different from its INS spectrum. This difference 
is schematically described in Fig. [5] and Fig. [7] 



VII. WEAK MAGNETIC FIELD 

Finally, let us briefly discuss the effect of weak mag- 
netic field. First we consider magnetic field H z applied 



c 

CD 

E 

(U 
X 
*i— 

4— I 

E 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 

a 

FIG. 5. The a-dependence of the matrix element 
Q (one t + ,one t~\V. R \0} a - (N = 16) The state 
|one t + , one t~) a corresponds to the scattering state of a soli- 
ton and an antisoliton in the vicinity of the critical point 
a = q c . Insets: Extrapolation of the matrix element at 
a = 0.45 to N — > CO. 







UJ 



FIG. 6. Schematic view of the INS peak positions of NTENP. 
The first peak corresponds to the soliton and the antisoliton, 
and the second peak does to the breather. 



in z direction. Total magnetization ^ . Sj is then a con- 
served quantity. Thus, the dispersion of each elementary 
excitation are simply shifted according to its magnetiza- 







FIG. 7. Schematic view of the RS peak positions of NTENP. 
The sharp peaks correspond to a single elementary excitation 
are absent due to the selection rule. 
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tion as 



e s = \l k 2 v s 2 + M s 2 - H z (soliton), (33) 

(.§ = \J k 2 v s 2 + M s 2 + H z (antisoliton), (34) 

ei = \Jk 2 v s 2 +Mi 2 (1st breather). (35) 

Phase transition occurs when the bottom of the soliton 
dispersion touches zero at H z = Ms- 

Next we consider the effect of magnetic field H x ap- 
plied perpendicular to the anisotropy axis z. In this 
case, the total magnetization ^ . 5J is no longer con- 
served. Therefore, a rearrangement of elementary exci- 
tations would occur. The mixing may be described by 
the effective 3x3 Hamiltonian in the one-soliton, one- 
breather, one-antisoliton subspace: 

(M s \ „ / V2 \ 
n cS = bM s - K/2 y/2\, 
\ M S J 1 \ y/2 / 

(36) 

where b is the mass ratio (Tl0|) . 

Diagonalizing this, we obtain the masses of 3 elemen- 
tary excitations in the transverse magnetic field as 

M A = M s , (37) 

(1 + b)M s ± JM s 2 (b-l) 2 + 4H x 2 
M± = y — • (38) 

Namely, one of the elementary excitations (which is an 
antisymmetric linear combination of soliton and antisoli- 
ton) has constant mass Ms- The other two masses de- 
pend nonlinearly on H x . The lowest mass M_ is reduced 
to zero at a critical field, where a phase transition occurs. 

This is essentially identical to the preceding analysis 
of anisotropic Haldane chains in Ref. l28| which was ap- 
plied to NDMAP 3 ? and NTENPi^, where the in-plane 
anisotropy E is also taken into account. As pointed out in 
Refs.[l5landli"7l the scattering intensity corresponding to 
the higher-energy elementary excitation (with the M + ) is 
suppressed in a moderate transverse field, where the ele- 
mentary excitation is absorbed into the excitation contin- 
uum. Systematic analysis of the effects of the transverse 
field in the SG theory framework is left for future. This 
will be necessary for discussion of ESR spectra^ in the 
present approach. 

VIII. SUMMARY AND DISCUSSION 

In this study, we found that the low-energy excita- 
tions of NTENP arc well described in the framework 
of the SG field theory. The elementary excitations of 
NTENP correspond to the soliton, the antisoliton and the 
breather, which are respectively created by applying 5+, 
5~ and Sq on the ground state. Their correspondence 



with the original spin model is deduced from the dimer 
limit a = 0, where the original model is exactly solv- 
able, and the numerical exact diagonalization calculation 
which shows that the excitation structure is smoothly 
connected from the dimer limit a = to a ~ a c . This 
picture of excitations well explains the preceding INS ex- 
periments on NTENP and numerical calculations of the 
DSF. 

Based on the established picture, we found that the 
sharp peaks, each of which corresponds to a single 
elementary excitation, vanish in the RS spectrum of 
NTENP. This is qualitatively distinct from the INS spec- 
trum. A slight disagreement with the experiments may 
be caused by the neglect of the E term anisotropy in 
the model Hamiltonian @, and the use of the effective 
field theory, which is only asymptotically exact in the 
low-energy limit, 

Although our primary focus in this paper was NTENP, 
the present result could be applied to other systems in 
the dimer phase of the Hamiltonian ([2]). Depending on 
the value of the uniaxial anisotropy D, the number of 
breathers changes: The system has two breathers for < 
D < 0.1J, one breather for 0.1J < D < 0.5J and no 
breather for 0.5 J < D. 

Besides the bond-alternating chain ([2]), there are also 
completely different class of systems that are described by 
the SG field theory. An example in quantum magnetism 
is the 5 = 1/2 antiferromagnetic chain with a staggered 
field. It describes several materials in a magnetic field, 
such as Cu benzoate. While this system is physically 
quite different from the model discussed in the present 
paper, they are both described by SG theory^ The dif- 
ference appears in the correspondence between the spin 
operators and the boson field. In fact, in the staggered 
field case, soliton and antisoliton couple to 5 Z while n- 
th breather couples to S x if n is even, and to S v if n is 
odd. (Here we define the direction of the staggered field 
as x direction.) As a consequence, observable spectra are 
quite different from the S=l chain studied in the present 
paper. 

As a final remark, in the context of the SG field the- 
ory, difference between the dimer phase and the Haldane 
phase is nothing but the sign of the cosine term in the 
Lagrangian density ([5]). This change should not yield 
a qualitative difference in the properties of the SG the- 
ory, although its physical consequences may be altered. 
The SG field theory could provide a starting point for a 
unified understanding of dynamics in different phases of 
5 = 1 antiferromagnetic chains. In a related but some- 
what different context, anisotropic 5 = 1 chains without 
bond alternation was studied in terms of SG theory^! It 
would be interesting to compare the two systems in more 
detail. 
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